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Dynamics of particle production by strong electric fields in non-Abelian plasmas 
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We develop methods for computing the dynamics of fermion pair production by strong color 
electric fields using the semi-classical Boltzmann-Vlasov equation. We present numerical results for 
a model with SU(2) symmetry in (1+1) Cartesian dimensions. 
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I. INTRODUCTION 

High energy collisions of heavy ions at RHIC, and soon 
to be carried out at the LHC, are thought to produce 
a quark-gluon plasma (QGP), yet efforts to model this 
plasma using quantum field theory methods have proved 
to be particularly difficult. Since heavy ion collisions 
are dynamic events, finding appropriate approximation 
schemes are important. The QCD-SU(3) case, for ex- 
ample, requires the use of three quark colors for each 
flavor and eight gluon fields. The non-Abelian charac- 
ter of QCD, preserving gauge invariance, and dynamic 
renormalization schemes all add to the difficulties. 

In a recent paper, we have demonstrated that clas- 
sical theories based on Boltzmann-Vlasov (B-V) trans- 
port equations have worked well for the QED prob- 
lem in (1+1) and (3+1) dimensions in both Cartesian 
and boost-invariant coordinates. So it seems reasonable 
to try to develop a transport equation for non-Abelian 
gauge theories. Fortunately, such transport equations 
have been worked out for some time, using the Wong 
equations pQ. (See, for example, Refs. [3] and [3| and ref- 
erences therein.) In this paper, we apply B-V transport 
equation methods to find the classical distribution func- 
tion for particles with SU(2) internal symmetry with an 
initial large electric field in the z-direction and no par- 
ticles present at t — 0. Particles are created by means 
of the Schwinger mechanism for particles with internal 
symmetry. We derive the necessary equations in (3+1) 
dimensions in Cartesian coordinates, but only give nu- 
merical results in (1+1) dimensions. The methods used 
here can be extended to boost-invariant coordinates in 
(3+1) dimensions and SU(3) symmetry. 

The Wong equations have also been used to develop a 
theory of non-Abelian fluid dynamics. See Refs. [3] and [51 



and to study unstable quark-gluon plasmas [BJ. 

Similar work for the SU(2) case using a Wigner func- 
tion has been carried out by Skokov and Levai [7] using 
methods described in Prozorkevich, et al. [5]. 

We discuss the QFT problem in Section [n] In Sec- 
tion Ini] we introduce the Wong equations for a particle 



with non-Abelian internal symmetry and find conserved 
quantities for the case we consider. In Section |IV[ we 
define particle distribution functions and in Section [V] 
we write down the B-V equation for our case. In Sec- 
tion [VlJ we give the Schwinger formulas for particle cre- 
ation, modified for our non-Abelian case, and in Sec- 
we find a solution of the B-V equation 



tion 



VII 



In 



Section VIII| we derive polarization currents. Numerical 
results are given in Section [Xj and conclusions in Sec- 
tion ED 



II. THE BACK-REACTION PROBLEM IN QCD 

The colored quarks obey the Dirac equation, 

{ ^ [tfy +gAft(x)] — M} tJj(x) = , (2.1) 

where A^(x) — A®(x)T a . The generators obey a Lie 
algebra: [T a ,T b ] = if abc T c . The sign of g here is 
chosen to correspond to that in Ref. [H In this paper 
we consider the situation in Cartesian coordinates where 
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<7 M!/ = diag(l,— 1, — 1,— 1). From Dirac's Eq. (2.1), the 
color current, 

J a >»(x)=g{j(x)^T a yb(x)), (2.2) 

obeys a conservation relation, D® ,c (x) J c '^(x) = 0, where 

D^{x)^5 ac d^+gf abc Al{x). (2.3) 

For the gauge field, the gauge field tensor F fiU (x) — 
F ^{x)T a is given by 

K,{*) = dMx)-d v Al{x)+gf abc A b (x)Al{x). (2.4) 
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and satisfies the equation of motion, 

D a ^ c (x)F c <^(x) = {J a ' v ) 



(2.5) 



In this last equation, we have taken the expectation value 
of the Dirac color current in some initial state of the sys- 
tem. The gauge field is treated as a classical commut- 
ing field. Eqs. pi} and pl>) constitute the QCD-QFT 



problem with feedback. 

In order to simplify this problem, we choose a gauge 
and require that the field A a ' fi (x) be in the z-direction 
and depend only on t for all a. That is, we put 



A a ^(x) = (0,0,0, A a {t)) . 



(2.6) 



Then, from (2.4), the only non- vanishing terms are given 
by 

F« z (t) = -F» t (t) = -d t A a (t) = E a (t) . (2.7) 

From ( |2.5[ ), for fj, = t and v — z, we find 

d t F a < tz {t) = -d t E a {t) = J a ' z {t) . (2.8) 

We can also have ji = z and v — t, in which case we find 

g f abc A b z (t) F c ' z \t) = g f abc A\t) E c (t) = J M (i) . (2.9) 

We can write these field equations in an obvious vector 
notation as 

d t A(t) = -E(t) , (2.10a) 
d t E(t) = -3 z (t), (2.10b) 
gA(t)xE(t)=3 t (t), (2.10c) 



Taking the derivative of p. 10c ) with respect to t and 
using (2.10b) gives a conservation equation: 



d t 3 t (t)+gA(t) x 3 z {t) =0. (2.11) 
The gauge field energy-momentum tensor is given by 

&"(x) = \g^F- p {x)F a ^{x)+F a ^{x)g a0 F a ^{x), 



and using the field equations p.5[ ), the divergence of the 
field tensor reads 



d fl e»»(x) = -F v "(x)-j a (x) 



(2.12) 



For the case we consider here, the field tensor can be 
written as 



E 2 (t) 



diag( 1,1, 1,-1). 



(2.13) 



The conservation equation (2.12) reduces to the two 
equations 

d t [E 2 {t)/2] = E(t) ■ [d t E(t)} = -E(t) • 3 z (t) , (2.14a) 



= E(t) ■ J*(t) , 



(2.14b) 



both of which are a result of dotting E(i) into the equa- 
tions of motion (2.10). 



III. CLASSICAL PARTICLE EQUATIONS 

A classical microscopic theory of colored charges in- 
teracting with a non-Abelian gauge field has been devel- 
oped by Wong in 1970 pQ. We follow the development in 
Litim and Manuel 0. The Wong theory for SU(N) as- 
signs to each particle on a trajectory x^^s) in space-time 
N 2 — 1 color charges Q a (s) which depend on the trajec- 
tory proper time. The color charges are to be thought 
of as replacing average values of the group generators, 
Q a (s) H> (T a (s)), in the Heisenberg representation. 

The proper time interval ds is defined by 



(ds) 2 = 9fiU (x) dx^dx" , 



(3.1) 



and the four velocity w M (s) and kinetic momentum /c M (s) 
along the path is defined by 



dx^(s) 

k"(s) = Mu tl (s) = M^-^ 

ds 

= {kt(s),k x (s),k y (s),k z (s)) . 



(3.2) 



Instead of starting with a Lagrangian, we simply state 
the Wong equations of motion for non-Abelian fields, 



, r dkn(s) 



M 



ds 
dQ a (s) 



gF^(s)k»(s)Q a (s), 
-gf abc A b Js)Q c (s)k^s) 



Writing Q a (s) as a function of x M (s), we have 

dQ^s) _ dQ?(x) 



u"(s) 



dx» 



(3.3a) 
(3.3b) 

(3.4) 



so that the equation of motion for Q a (s), Eq. (3.3b), can 
be written in terms of the covariant derivative, 



k»D™{x) Q c (x) = 0. 



(3.5) 



where D® c (x) is given in (2.3). 

In a vector notation, Eqs. (3.3) can be written as 



M 



M 



dfc M (s) 

ds 
dQ(,s) 



= gQ(s)-F^(s) k v (s), 



= - 9 A,( s )xQ( s )P( S ), 



and (3.5) becomes 



k^D^Qix) = 



(3.6a) 
(3.6b) 

(3.7) 



where D^(x) is the covariant vector derivative operator 
D jI (i) = a (I + s A fI (.i)x . (3.8) 



Eq. (3.6b) says that the color charge vector Q(s) pro- 
cesses about the vector A u ( s) u^j s) along the trajectory 
in space-time. Multiplying (3.6a) by fc M (s), and requir- 



ing the momentum to be on the mass shell, leads to the 
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fact that the length of the kinetic momentum vector is 
conserved, and given by 

^(s)^(s) = k 2 (s)-k 2 x (s)-k 2 Js)-k 2 z (s) =M\ (3.9) 



IV. DISTRIBUTION FUNCTION 

We define a particle distribution function f(x, k, Q) 
such that the average convective particle color current 
density J£ on (s) is given by 



k t (s) = ^kl(s) + kl(s) + M 2 , (3.10) ^ m (x)=gjDkDQk^Qf(x,p,Q) 



(4.1) 



where kj_(s) — k 2 (s) + k 2 (s). Dotting Q(s) into 
Eq. (3.6b I leads to the observation that the length Q 2 of 



the color charge vector, which is the quadratic Casimir, 
is conserved. For SU(3), the cubic Casimir is also con- 
served. 



For our case, Eq. (3.6a) becomes 
dk t (s) 



M 
M 



ds 
dk z (s) 
ds 



gQ(s) -E(s) k z (s) , 
gQ{s) -E(s) k t {s) , 



(3.11a) 
(3.11b) 



with k x and k y (and consequently k±) constants of the 
motion. For our case, Eq. (3.6b I becomes 



M dQ( f ) ^ 
ds 



The dot product of A(s) with Eq. (3.12) gives 



MA( S ).¥ = ; 



and since 



/ x / x dt dA(t) dA(s) 
fc t (s)E(s) = -A/- ^ = -Ai- 



ds dt 



ds 



Eq. (3.11b) can be written as a total derivative: 



M-[k z (s)+gQ(s)-A(s)]=Q. 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



So the quantity in brackets is a constant of the motion, 
which we call p z . Then we have 



by 



k z (s)=p z -gQ(s)-A{s). (3.16) 
For a single particle, the color current J^(x) is given 



J"(x)=g J dsQ(s) 



^Ms^x-xis)}, (3.17) 



and the matter energy momentum tensor t^ v (x) by 
tT(x) = / ds 



dx " ( - S h ,J (s)S 4 [x~x(s)}. (3.18) 



ds 



In the next section, we introduce a distribution function 
for an ensemble of particles. 



Similarly, the particle energy-momentum tensor density 
is given by 



(4.2) 



tf a/ (x) = J DftDQA"*" f{x,p,Q) . 

Here the momentum measure Dfc is given by 
2re{k t )S(k 2 - M 2 )d i k 



Dfc = 



(2tt) e 



(4.3) 



where r is a degeneracy factor which counts the number 
of species. For one flavor of fermions and anti-fermions 
with no spin r = 2, so in (1+1) Cartesian dimensions 



Dfc = 



dk z 



(4.4) 



In (3+1) dimensions for one flavor of fermions and anti- 
fermions with spin r = 4, and this factor is given by 



Dfc = 



k±dk± dk z 

7T 2 Uk±,k* 



"k ± ,k z = Jkl + kl + NP . (4.5) 



In subsequent sections of this paper, we work out the nec- 
essary equations in (3+1) dimensions — the translation 
to (1+1) dimensions essentially means that we omit the 
k± integration, set fcj^ = 0, and multiply the currents by 
7r. For SU(2), the color measure is 

VQ = c R d 3 Q5(Q 2 -q 2 ) (4.6) 

where the delta-function expresses conservation of the 
quadratic Casimir. Here cr is a normalization factor, cr 
and qi are set by the conditions, 



(4.7) 



Tr[l] =2 i — ^ / DQ = 2 , 
Tr[T-T]=3/2 i-> JbQQ-Q = 3/2 

Writing Q in spherical coordinates, 

Q 1 = J sin0 cos</> , 
Q 2 = J sin 9 sin cj> , 
Q 3 = J cos9, 



(4.8) 



we find from Eqs. ( |4.7[ ), the results: cr = 2/(ttv3), q% = 
3/4, and J = v3/2. Similar results are found for SU(3) 
where there are two conserved Casimir invariants [2J. 
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V. BOLTZMANN-VLASOV EQUATION 

The Boltzmann-Vlasov equation can be derived in 
parametric form by considering the total derivative of 
a function f[x(s), k(s), Q(s) ]. Using the equations of 
motion (3.3), we find 



M 



df[x(s),k(s),Q(s)[ 



M 



dx^(s) d 
ds dx^ 



dfc M (s) d dQ a (s) d 



dk., 



dQ a 



f(x,k,Q) 



(5.1) 



= W(s)B li [A](s)f(x,k,Q), 

where the B-V differential operator S[A](s) is defined by 

Bfi[A](s) = V^A] - g Q ■ F M „ (x) d kv , (5.2) 
V^A] = d„ -gA„(x) • Q x d Q . (5.3) 

Here 2^ [A] is a color-covariant derivative operator, in- 
variant under gauge transformations. So now requiring 
that k to be on the mass shell and Q satisfy the Casimir 
relations, the Boltzmann-Vlasov equation is given by 



A" B„[A]{s) f(x, k, Q) = fc* C{x, k, Q) , 



(5.4) 



where C(x,k,Q) is a source term. For our case, the dis- 
tribution function is a function of t, k z , and Q only. 
Then the Boltzmann-Vlasov equation (5.4 1 becomes 



d_ 
dt 



f abc A b (t) Q c 



d 



(5.5) 



gQ a E a (t) 



d 
dk~ z 



f(t,k ± ,k z ,Q) = C(t,k ± ,k z ,Q) 



A. Conservation of color current 

We first study the divergences of the convective color 
current and the particle energy-momentum tensor. For 
the convective current, we multiply (5.4 1 by gQ and in- 
tegrate over Dk and DQ. This gives 

gjvkj DQ W Q B^ [A] f{x, k, Q) (5.6) 

= g JmJ DQQk*C(x,k,Q). 



In the first term on the left-hand side of (5.6), the deriva- 



tive with respect to x M can be factored out of the integral. 
The second term gives a factor g A^(x) xJ£ on (x) by parts 
integration, and the third term vanishes by parts integra- 
tion and the antisymmetry of F M „, yielding the equation 



n^(x)J^ n (x)=gV(x), (5.7a) 
V(x) = jDkDQu kxikm QC(x,k,Q). (5.7b) 



where D M (i) is the covariant vector derivative opera- 
tor, defined in Eq. ( |3.8| . For our case, the convec- 
tive color current four- vector is of the form J^ on (t) — 
(j£ on (0,0,0,J c z on (i)), where 



,(*) = 



o r~\-oo p 

k ± dk ± dk z DQ 

J — oo J 

x Qf(t,kj_,k z ,Q) , 

o p-\-oo p 

k±dkj_ dk z DQ 

J — oo J 

■Qf{t,k±,k*,Q) 



(5.8a) 



(5.8b) 



u} k,,k z 



So for our case, Eq. (5.7a) becomes 

d t 3 t con (t)+gA(x)x3 z con (t)=gV(t) 



(5.9) 



The complete current from the gauge field equations must 
satisfy Eq. (2.11). We shall see below how to define a 



polarization current which, when added to the convec- 
tive curre nt, will produce a total current which satisfies 
Eq. (|2.11). 



B. Conservation of energy and momentum 

To find the conservation law for the particle energy- 
momentum tensor, we multiply ( 5.4 ) by k v and integrate. 
This gives 

J Dk J DQk^k v B^.[A]f{x,k,Q) (5.10) 

= J Dk J DQ k v k* C(x, k, Q) . 

Again the derivative in the first term on the left-hand 
side comes out of the integral. By parts integration, the 
second term vanishes but the third term yields — g Q • 
F va {x) k a , so we find 

0„ t» v (x) = F va {x) ■ J c ° n {x) + G u {x) , (5.11a) 

G"(x)= [[dk] ( DQ k v C(x, k, Q) . (5.11b) 



Following a method introduced by Gatoff, Kerman, and 
Matsui [9] , we will find in Section VIII below that for our 
case, we can write 



G u {x) =F' MJ {x) • JP ol (x) , 



(5.12) 



where J? 01 (a:) is a polarization current. Defining the total 
particle current as J (r (x) = J™ n (x) + JP ol (x), Eq. pTTT] ) 
becomes 



d^{x) = F™(x)-3 a {x). 



(5.13) 



The field energy-momentum tensor satisfies Eq. ( |2.12[ ), 
where J CT (x) is the total current, so the divergence of the 
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sum of the matter and field energy momentum tensors, 
T^(x) = t^(x) + e» u (x) is given by 



d fi T^(x)=0, 



(5.14) 



and is conserved. For our case, the particle energy- 
momentum tensor ^"(x) is diagonal, and given by 

t^(t) = d^g(e(t),p x (t),p y (t), Pz (t)) , (5.15) 

where the particle energy and pressure densities are given 
by 



e(t) 



fJ 



oo /*+oo /> 

k ± dk ± dk z DQ 

J —oo J 

x ^k ± ,k z f(t,k±,k z ,Q) 

oo 



I) 



Pz(t) = 



±dfc_L J 


dk z / DQ 




OO J 


K ± 


f(t,k±,k z ,Q) 


UJ k±,kz 


j_dk ± J 


dk z / DQ 




OO J 


kl 


f{t,k±,k z ,Q) 


UJ k±.k z 



(5.16a) 



(5.16b) 



(5.16c) 



with p x (t) = p y (t) and p ± (t) = p x (t) + p y (t) . 

VI. PARTICLE CREATION 

A. QFT pair production 

Pair production rates via the Schwinger mechanism by 
strong non-Abelian gauge fields have recently been calcu- 
lated using a one-loop approximation in QFT for bosons 
by Nayak and Nieuwenhuizen [10] and for fermions by 
Nayak [IT] . The key to the calculation was to diagonalize 
the matrix E = J2 a E a T a . Let us write the eigenvalue 
equation for the matrix E as 

Ee i = E i e i . (6.1) 

Then the rate of fermion pair production in (3+1) di- 
mensions is given by a sum over the eigenvalues |llj . 



dK 



d 7 N 



dtd 3 xd 3 k 

\gE, 



7r(fcj + M 2 ) 

\gE~\ 



(6.2) 



where M is the quark mass. For the case of SU(2), the 
eigenvalues are given by E± — ±|E|/2, so the two eigen- 
values give identical contributions. 
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FIG. 1: Plot of A(t) as a function of t for a solution of the 
B-V equation with feedback for SU(2) gauge fields. At t — 0, 
we took A(0) = and E(0) = (4, 2, 1), with M = 1 and g = 1. 
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FIG. 2: Plot of E(t) as a function of t for a solution of the 
B-V equation with feedback for SU(2) gauge fields. Same 
initial conditions. 



B. Classical pair creation rate 

For the classical calculation, it would seem natural to 
map E a (t)T a H> E(t) • Q, and replace the sum over eigen- 
values by an integral over Q. This mapping produces a in- 
dependent particle production rate C(t, kj_, k z , Q) given 
by 

C(t,k ± ,k z ,Q) = \gQ-B(t)\R(t,k^Q)S(k z ) , (6.3) 
R(t, k±, Q) = Pit, k±,Q) S(t, k x ,Q) , 

where 

(6.4a) 



P(t,k±,Q) = l-2f (t,k±,Q), 

7T (fcj + M 2 ) 

' lffQ-E(t)| 



S(t,k±,Q) = -{ 1 -exp 



.4b) 



and where we have set f (t, k±,Q) = f(t, k±, 0, Q). Here 
we have multiplied (6.2) by S(k z ), since most of the par- 
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FIG. 3: Plot of Jcon(i) as a function of t for a solution of 
the B-V equation with feedback for SU(2) gauge fields. Same 
initial conditions. 
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FIG. 4: Plot of Jp i(i) as a function of t for a solution of 
the B-V equation with feedback for 5*17(2) gauge fields. Same 
initial conditions. 



Here s represents the path-length in space-time and x(s'), 
k(s') and Q(s') are solutions of the particle trajectory 
equations ( |3. 11 1 for values of s between s and s. The 
trajectory equations must be integrated backwards, start- 
ing with the "current" values of x, k, and Q, and winding 
up with a value of sq such that i(so) = 0. We arbitrarily 
choose s = so the values of s' are negative. Then the 
"initial" conditions for the trajectory solutions are 



x(0) = x , fc(0) 



0(0) 



(7.2) 



which are the current values of x, k, and Q, and the 
final value of So defined by t(so) = 0. Since no particles 
are present at t = 0, /[0, fe(so), Q(so) ] = 0, and (7.1) 
becomes 



f(t,k_L,k x ,Q) 



(7.3) 



= / ds' k t (s')C[s' ,kx,k z (s'),Q(s')}/M 



and using (6.3), we find 



f(t,k±,k z ,Q) = f ds'k t (s') (7.4) 
x | gQ(s') ■ E(a') | R[ s', k ± , Q(s') } 6[ k z {s') }/M . 



From (3.16), we have 



k z (s')= Pz -gQ(s')-A(s') (7.5) 
= k z + g[Q-A(s)-Q(s')-A(s')]. 



Inserting this into the delta-function in ( 7.4 ) , and noting 
that 



M 

we find that 



dk z (s') 
ds' 



, 9 Q( S ')-E( S ')fc t (s'), 



(7.6) 



f{t,k ± ,k z ,Q) 



(7.7) 



n,k±,Q(s n )]&[t(s n )]G[t-t(s n )[ 



ticle production occurs at k z = 0, and introduced a Pauli 
supression factor at k z =0. 



VII. SOLUTION OF THE 
BOLTZMANN-VLASOV EQUATION 

In this section, we solve the BV equation using the 
method of characteristics. From Eqs. (5.1) and (5.4), a 



solution of the BV equation can be found by integrating 
the source term over a classical particle path trajectory 
from so to s, 

f[x(s),k(s),Q(s)]=f[x(s ),k(s ),Q(s )] (7.1) 
+ f ds' k t (s')C[x(s'),k(s'),Q(s')]/M , 

J s 



where s n is a solution of the equation 

k z + g[Q-A(s = 0)-Q(s n )-A(s n )] = 0, (7.8) 

and s n must be in the range Sq < s n < 0. This completes 
the solution using the method of characteristics. 

Since the Pauli term depends on the distribution func- 
tion evaluated at k z — 0, Eq. (7.7) needs to be solved 
explicitly for the case when k z — 0. For this case s n is a 
solution of the equation 



Q • A(s = 0) = Q(s n ) • A(s„) 



(7.9) 



One such solution will always be s n = s = 0, or t n = t. 
Then using the relation 9(0) = 1/2, and solving (7.7) for 
f (t,k±,Q) = f(t,0,k±,Q) gives 



fo(t,k±,Q] 



S(t,k±,Q)/2 + Z(t,k±,Q) 
l + S(t,k±,Q) 



(7.10) 
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where 

Z{t, k±,Q) = R ( s ^ k ^Q) ■ ( 7J1 ) 

s n <0 

This completes the solution of f(t,k±,k z ,Q) using the 
method of characteristics. The method used here for 
non-Abelian symmetries differs from that used for QED 
since in this case we must numerically solve the trajec- 
tory equations for Q(s) given the final values of x, k. 
Here the trajectory solution for Q(s) does not have a 
conserved quantity as we found for k z (s) in Eq. (3.16|. 
Note that we do not need to have the values of k or Q 
at s — so since by the backwards integration, we find the 
"initial" values needed to obtain the final values at t. 

However, it is not necessary to find the complete dis- 
tribution function f(t,k±,k z ,Q) in order to compute 
the currents since we can make use of the (5-function in 
C{t, k±, k z ,Q) to replace the integral over k z by an in- 
tegral over s' . This is a great advantage since it means 
that we only need to find the special distribution func- 



tion fo(t,k±,Q). Substituting (7.4) into Eqs. (5.8) and 



integrating over k z , the convective color currents can be 
found from the equations, 



^ J k ± dk ± J DQ Q J ds' kt(s') (7.12a) 



\gQ( 3 ')-E( S ')\R[ S ',k ± ,Q(s')]/M., 



,(«) = 



poo p rO 

% / kxdkx / DQ Q / ds' kt(s') (7.12b) 

x k *(°'> s '® | g Q( S ') • E( S ') | R[ s>, k x ,Q(s') }/M , 

where we have put 

k z (s', s,Q)=g [ Q(s') • A(s') - Q • A(s) ] (7.13a) 



u k± {s', s, Q) = ^kl + k 2 z (s', s ,Q) + AP . (7.13b) 



For the energy and pressures, substituting (7.4) into 
Eqs. (5.16) gives 



e(t) = 



^ jf k±dk± J DQ^ds'fc t (s') (7.14a) 



x Lo k± (s', s,Q)\g Q(s') ■ E(s') | R[ s', k x ,Q(s') }/M , 

POO /> nO 

PJ-W = 4 / fc ^ d ^ / D Q / ds ' fc t( s ') (7.14b) 



k\ 



w fci (s',s,Q) 

„ fc2( s ', s ,g) 



| 5 Q( S ')-E( S ')|i?[ S ',^,g( S ')]/M, 



J^jf k ± dk ± J DQ ds'k t (s') (7.14c) 
, lfl Q(a , )-E( fl / )|iJ[ a / J * X| Q( a , )]/M. 

w fci (s',s,g) 

In Eqs. ( |7.12[ ) and (7.14), we only need to find the dis- 
tribution function fo(t, k± , Q) to obtain the currents, en- 
ergy, and pressures. However we will still need to back- 
integrate the trajectory equations to find Q(s')- 




14 16 18 20 



FIG. 5: Plot of the particle, field, and total energies as a 
function of t for a solution of the B-V equation with feedback 
for SU(2) gauge fields. Same initial conditions. 
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FIG. 6: Plot of the particle production rate dn/dt and total 
particle production n(t) from Eq. (9.1 1 as a function of t for a 
solution of the B-V equation with feedback for SU(2) gauge 
fields. Same initial conditions. 



VIII. POLARIZATION CURRENTS 

The total current is the sum of the convective and po- 
larization currents, 



n*)=J£on(*) + J£ol(*)> 



and satisfies the gauge field equation Eq. (2.11), 
atJ*(t)+ffA(i) xJ z (t) = 0. 



On the other ha nd, t he con vect ive part satis fies Eq. ( |5.9[ ) . 
So substituting (8.1) into (8.2) and using (5.9), we find 



(8.1) 



(8.2) 



that the polarization current satisfies 

d t Jl oi (t) + gA(t) x j; ol (t) = -gV(t) , (8.3) 



which fixes J pol (t) in terms of Jp 0l (£) and V(i). Substi- 
tuting (6.3 1 into Eq. (5.7b), V(i) is given by 



V W = / fc ^ d ^ / D( 9 Q 1 .9 Q • E(t) I fl(t, fc x , Q) • 
The current component Jp Q i(i) is obtained from energy 



conservation. Substituting (|6.3| into Eq. (5.11b I gives 



G\t) = Ay fc^J DQ (8.5) 

xw fcxi0 |Q-E(i)|i2(t,Ax,Q). (8.6) 
Now we can always write, 

|Q-E(i)|=sga[Q.E(t)]Q.E(t). (8.7) 
So G*(t) = E(t) • Jljt), where 



JSoi(i) 



/c_l dfc ± / DQ 



x Lj k±t0 sgn[ Q • E(i) ] Q JZ(t, fe±,Q) 



So the rate of change of particle energy density is given 
by 



(8.9) 



where 3 z (t) is now the total z-component of the current. 
Adding the field energy density from Eq. (2.14a I to the 



particle density, we find that the total energy density is 
conserved: 



d t [e(t) + E 2 (t)/2] =0. 



(8.10) 



So in this section, we have found equations for the po- 
larization contribution to the current by requiring that 
total energy, including field energy, be conserved for our 
source terms. The polarization current is a required mod- 
ification of the gauge field equations to account for the 
creation of particle and anti-particle pairs. 
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FIG. 7: Plot of the average value of Jcon(i) as a function of 
t for a solution of the B-V equation with feedback for SU(2) 
gauge fields. Same initial conditions. 



0.0010 




-0.0010 



10 12 14 16 18 20 
t 



FIG. 8: Plot of the average value of Jp 0l (i) as a function of 
t for a solution of the B-V equation with feedback for SU(2) 
gauge fields. Same initial conditions. 



NUMERICAL METHODS 



IX. PARTICLE PRODUCTION 



Integrating the B-V equation (5.4) over k z and Q gives 
dn(t, k±) 



Of 



J^-\gQ--E(t)\R(t,k^Q), (9.1) 



where the particle density n(t, k±) is given by 
"+°° dk z 



n(t,kj_) = 



r-t-°° AU f 

J J VQf(t,k ± ,k Zl Q). (9.2) 



A plot of the particle density n(t, 0) as a function of time 
for the case of (1+1) dimensions is shown in Fig. [6] 



Let us review the equations we need to solve, 
fields satisfy the equations 



The 



d t A(t) 



-E(i) , 

-r{t) , 



(10.1a) 
(10.1b) 



where the current is given as the sum of convective and 
polarization currents, J M (i) = JconW + Jpoift)' Compo- 



nents of the currents are given by Eqs. (7.12) and (8.8) 



Jpol(i) is given by the solution of Eq. (8.2), which is 



stepped out along with Eqs. (10.1). Back integrations 



of the trajectory equations of motion are required to find 
kt(s'), k z (s'), and Q(s') for the currents. The trajectory 
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equations are given by 
dkt(s) 



M 
M 
M 



ds 
dk z (s) 

ds 
dQ(s) 

ds 



gQ(s) ■ E(s) k z (s) , (10.2a) 
gQ(s) ■ E(s) k t (s) , (10.2b) 
g A(s) x Q(s) k z (s) , (10.2c) 



with 



M ^) = fct(s ) , M ^ = jfe,(a) . (10.3) 



ds 



ds 



The backward integration is started at a proper time 
point s = where t(0) = t. At this point, particles are 
created with zero z-component of momentum, k z {t) = 0. 
So given values for fcj_(0) = k± and Q(0) = Q, and using 
the parametric equations, the "initial" condition at s = 
for the backward integration is then specified as 



t(0)=t, k t (0)-- 
z(0)=0, k z (0)=0 



k\ 



M 2 



(10.4) 



Eqs. ( 10.2 ) and ( 10.3 1 are integrated backward to t = 0, 
using the known values of A(s') and E(s') for Sq < s' < 0, 
obtaining and values for fct(s'), k z (s'), and Q(s'). Using 
these values, the integral over s' for the currents, energy, 
and pressures can be done, and the process repeated for 
all values of k± and Q needed for the integral over s'. 
During this back integration, it is also necessary to find 
the special distribution function fo(t, k±, Q) at each step. 
This is given by Eqs. (7.10) and \7.ll I and stored glob- 
ally. 

A standard fourth order Runge-Kutta routine is used 
to step out the field equations (10.1) and the back inte- 
grations (10.2). In order to improve the accuracy of the 



Runge-Kutta method for the field equations, it is desir- 
able to compute the derivative of the z-component of the 
current at each step, and use this for a linear interpola- 
tion for the currents during the steps. Since the polar- 
ization current is generally quite small, we use only the 
derivative of the z-component of the convective current. 
This is given by 

d t 3 z con (t) (10.5) 

= ^ k ± dk ± / DQ(fci + M 2 )Q[ 5 Q-E(i)] 

r° ,k t (s>) \ Q( S ').E( S ')\R[s',k±,Q(s')] 



ds' 



M 



ul ± (s',s,Q) 



The color particle number density is obtained from 
J*(f) = J* on (i) + Jp ol (t). The ^-component of the con- 



vective current is obtained from Eq. (7.12a) and the t- 



component of the polarization current is obtained by the 
solution of Eq. 



■ 2) 



Some results are shown in Figs. [T}|6] for the case where 
we have set M = 1, g = 1, and taken A(0) = and 



E(0) = (4, 2, 1) and ignored the Pauli correction. Here 
we took dt = 0.002 and ds = 0.0005 with 320 values of 
Q distributed in more or less equal area triangles about 
the surface of the Q-sphere. This seems to be sufficient 
for values out to t = 20. The program took about 24 cpu 
hours on a desktop-type machine. 

In Figs. [I] and [2] we see that the fields A(i) and E(t) 
oscillate with a period of approximately 20 in units of the 
mass to. A(i) remains negative during this time inter- 
val. The z-component of the convective current Jj: on (i), 
shown in Fig. [3j appears to saturate even though there is 
no Paui correction factor here. The polarization current, 
shown in Fig. [4j is about one-tenth the size of the convec- 
tive current but is necessary to create plasma oscillations 
and is required for energy conservation. In Fig. [5j we plot 
the particle and field energies. Energy is exchanged be- 
tween these two energies with the total energy conserved 
to an accuracy of about 1%. The rate of particle pro- 
duction and total particle production is shown in Fig. [6] 
Over half the total production is during the first few os- 
cillations of the fields. 

The time-components of the convective current should 
vanish since the integrand in Eq. (7.12a) is odd under 
Q 4 -Q. As a check, we show by direct calculation in 
Fig. [7] that this component of current is less than 5 x 
10~ 3 in magnitude for t < 20. The calculation begins 
to show instability for t > 14. Likewise V(i) vanishes 
since the integrand in Eq. (8.4) is odd under reversal of 



Q. Eq. (8.3) then reduces to 



d t J t Ut)+gA(t)xJ*(t) = Q. 



(10.6) 



Numerical solutions of ( 10.6 ) are shown in Fig. [8j Here 
we see that all components of Jp D i(i) are less that 1 x 10~ 3 
which, to the accuracy of this calculation, is consistent 
with zero. So the auxiliary Eq. ( |2.11 1 is satisfied by this 
calculation. 



XI. CONCLUSIONS 

In this paper, we have derived semi-classical transport 
equations for particle production by color gauge fields in 
a quark plasma with SU(2) gauge symmetry. We solved 
a particularly simple Cartesian (l+l)-dimensional case, 
where the electric field is constrained to be in the z- 
direction, and showed how to obtain numerical solutions 
to the B-V equation with a Schwinger-type source term, 
starting with no particles and an initial large color field. 
Our results conserve the average color and energy of the 
plasma, and satisfy the gauge field conditions. 

With additional computing power, our calculation can 
be extended to boost-invariant (3+l)-dimensions for the 
SU(2) case. It would also seem possible to extend the 
calculation to the more interesting SU(3) case, where we 
could explore particle production as a function of the two 
conserved Casimir invariants. 



10 



Acknowledgments 

This work was performed in part under the auspices of 
the United States Department of Energy. The authors 
would like to thank the Santa Fe Institute for its hos- 
pitality during the completion of this work. One of us 



(JFD) would like to thank Robert Carrier of the Univer- 
sity of New Hampshire Research Computing Center for 
helpful discussions. He would also like to thank Matthew 
Minuti for computing advice and Silas Beane for use of 
computing resources. 



[1] S. K. Wong, Nuovo Cim. 65A, 689 (1970). 

[2] D. F. Litim and C. Manuel, Phys. Rept. 364, 451 (2002), 

arXive:hep-ph/0110104v2. 
[3] B. P. Schenke, Ph.D. thesis, Johann Wolfgang Goethe- 

Universitat, Frankfurt am Main (2008), arXiv:hep- 

ph/0810.4306vl. 
[4] B. Bistrovic, R. Jackiw, H. Li, V. P. Nair, and S.-Y. Pi, 

Phys. Rev. D 67, 025013 (2003). 
[5] R. Jackiw (2003), arXive:hep-th/0305027v2. 
[6] C. Manuel and S. Mrowczyhski, Phys. Rev. D 74, 105003 

(2006). 



[7] V. V. Skokov and P. Levai, Phys. Rev. D 78, 054004 

(2008), arXive:0710.0229v2 [hep-ph]. 
[8] A. Prozorkcvich, S. Smolyansky, V. Skokov, and 

E. Zabrodin, Phys. Lett. B 583, 103 (2004). 
[9] G. Gatoff, A. K. Kerman, and T. Matsui, Phys. Rev. D 

36, 114 (1987). 

[10] G. C. Nayak and P. van Nieuwenhuizen, Phys. Rev. D 

71, 125001 (2005). 
[11] G. C. Nayak, Phys. Rev. D 72, 125010 (2005). 



